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PURPOSE:  This  Coastal  and  Hydraulics  Engineering  Technical  Note  (CHETN)  provides 
improved  predictive  equations  for  coastal  rubble-mound  stone-armor  stability.  The  engineering 
methods  outlined  apply  to  breakwaters  and  revetments  exposed  to  wave  loading.  The  stability 
equations  are  based  on  the  concept  that  the  maximum  wave  force  causing  armor  instability  is 
proportional  to  the  maximum  wave  momentum  flux  near  the  structure  toe.  This  concept  introduces  a 
more  physics-based  first  principles  approach  to  estimation  of  annor  stability.  The  new  equations 
explicitly  include  the  effects  of  nearshore  wave  height,  wave  period,  water  depth,  and  stonn  duration 
as  well  as  the  characteristics  of  wave  breaking  on  the  structure. 

INTRODUCTION:  Annor  stability  on  large  coastal  rubble  mounds  has  historically  been  based  on 
empirical  equations  that  relate  armor  buoyant  weight  to  the  maximum  wave  drag  force.  The  fluid 
velocity  in  the  maximum  drag  force  was  computed  using  the  shallow-water  wave  celerity 

C  =  sjgH  where  g  is  the  acceleration  of  gravity  and  H  is  the  wave  height  near  the  structure.  Hudson 
(1959)  and  other  early  researchers  used  this  formulation  to  develop  the  stability  number 
Ns  =  II/ADnso  as  the  basic  measuring  stick  for  stability.  Here  A  =  Sr  —  1  is  the  buoyant  density  of  the 
stone,  Sr is  the  stone  specific  gravity,  and  Dn 50  =  (Vnso)  where  V„5o  is  the  median  stone  volume. 
Van  der  Meer  (1987)  also  used  this  fonnulation  of  the  stability  parameter.  The  weakness  in  these 
methods  is  that  the  fluid  force  is  too  simplistic  to  be  generally  applicable.  For  shallow  water  depth- 
limited  applications,  the  maximum  wave  force  depends  on  the  local  water  depth  and  wave  period. 
Therefore  these  variables  should  be  included  in  the  derivation  of  the  equations. 

MOMENTUM  FLUX:  Melby  and  Hughes  (2004)  derived  a  formulation  for  armor  stability  by 
assuming  the  maximum  wave  force  was  proportional  to  the  maximum  wave  momentum  flux. 
Hughes  (2003a,  2003b,  2004,  2005)  introduced  the  wave  momentum  flux  for  coastal  applications. 
Assuming  irrotational  potential  flow  on  a  locally  flat  bottom  in  water  depth  h,  the  wave-averaged 
and  depth-integrated  radiation  stress  is  given  by  the  integration  of  wave  momentum  flux  over  the 
wavelength,  i.e., 


M  =  S„ 


L 


L  fix 

f  f(Pd  +P  wu2)dzdx 

0  -h 


(1) 


where L  is  the  wavelength,  r\x  is  the  free  surface  location,/;,/  is  the  dynamic  pressure,  pw  is  the  fluid 
density,  u  is  the  velocity  in  the  x-direction,  x  is  the  horizontal  coordinate,  and  z  is  the  vertical 
coordinate.  The  maximum  depth-integrated  wave  momentum  flux  is  given  at  the  wave  crest  as: 
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HI  2 

(Mf) max  =  /  (Pd  +  ?u2)dz 

—h 


Using  linear  wave  theory  values  for  u  and  pd  yields 


P  wgh2 


1  H  tanh  kh 

2  h  kh 
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(3) 


where  k  =  2i i/L  is  the  wave  number.  In  Equation  3,  the  first  term  on  the  right-hand  side  is  the 
dynamic  pressure  term  while  the  second  is  the  velocity  term.  In  general,  the  pressure  term  will 
dominate.  For  example,  for  low  steepness  waves,  the  velocity  tenn  will  only  contribute  5  percent  to 
the  maximum  momentum  flux.  For  waves  in  shallow  water  at  the  steepness  limit,  the  velocity  tenn 
will  provide  the  maximum  contribution,  roughly  30  percent  of  the  momentum  flux.  Equation  3 
assumes  waves  to  be  periodic,  linear,  and  sinusoidal,  and  it  ignores  the  momentum  flux  above  the 
still-water  level.  However,  in  shallow  water,  waves  are  nonlinear  with  peaked  crests  and  shallow 
troughs.  The  wave  forces  from  these  nonlinear  waves  can  be  very  different  from  those  estimated 
from  linear  waves.  Equation  3  will  underpredict  the  momentum  flux  under  nonlinear  wave  crests. 
For  MF/pwgh  =  0.2,  Equation  3  will  produce  an  underprediction  error  of  roughly  10  percent  while 
for  M/.jpwgh  =  0.8,  the  error  will  be  roughly  100  percent. 


The  maximum  wave  momentum  flux  increases  rapidly  for  nonlinear  waves  -  steep  waves  in 
shallow  water.  This  corresponds  to  the  case  where  armor  stability  is  at  its  minimum.  It  is  desirable 
to  develop  a  relation  that  can  characterize  the  stability  over  the  full  range  of  water  depths  expected. 
Hughes  (2003a,  2003b,  2004, 2005)  estimated  the  nonlinear  wave  momentum  flux  using  a  numerical 
Fourier  solution.  The  results  were  found  to  be  well  represented  by  the  empirical  equation: 


Mf 

,P  wSh\ 


An 


gT2J 


(4) 


with 


A0  =  0.639 


H 


2.026 


A  =0.180 


,  >*-0.391 


(5) 


Use  of  the  nonlinear  approximation  in  Equation  4  is  important  because  stability  is  at  its  minimum 
when  the  incident  wave  is  the  most  nonlinear.  The  actual  wave  momentum  flux  force  that  a 
particular  armor  unit  is  exposed  to  will  vary  from  the  value  given  by  Equation  4  due  to  the  slight 
error  in  the  numerical  approximation,  the  effects  of  the  sloping  foreshore,  armor  unit  position  in  the 
water  column,  and  armor  unit  extent  of  the  water  column.  Therefore,  a  stability  relation  derived 
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using  Equations  3  or  4  will  be  empirical.  Herein,  simply  assume  that  the  force  on  an  armor  unit  is 
proportional  to  the  depth-integrated  maximum  wave  momentum  flux  given  by  Equation  4  computed 
at  the  structure  toe. 

NEW  ARMOR  STABILITY  EQUATION:  Generalized  empirical  stability  equations  were 
developed  that  characterize  annor  instability.  These  stability  equations  were  fit  to  small-scale 
laboratory  data  from  Van  der  Meer  (1987).  The  data  correspond  to  limited  conditions  as  follows: 
mostly  nonbreaking  waves,  nonnally-incident  waves,  nearshore  slope  of  1V:20H,  mostly  non- 
overtopped  structures,  narrow  and  wide  armor  gradations,  permeable  and  impermeable  structures, 
homogeneous  structures,  and  angular  randomly-placed  armor  stone.  A  small  number  of  tests 
included  shallow-water  breaking  waves  and  an  additional  small  number  included  overtopped 
structures. 

Two  new  equations  resulted  from  the  fit  to  data  representing  the  two  most  important  breaker  types. 
Plunging  waves 

Nm  =  5.0(5 /  JV°5)°'2P°  l8Vcot0  s,  >  V  (6) 


Surging  waves 


Nm  =5.0(S/N^yzPul\cotQ) 


\0.5 -P 


-P/3 


.S'  <  ,S„ 


(V) 


where 


s  =  -0.0035  cot0  +  0.028 


(8) 


and 


N-~\  I  PS  <9) 

with  Ka=  1 .  Here  smc  is  the  critical  wave  steepness  on  the  structure,  P  =  structure  penneability, 
S  =  AJD1  nso  =  normalized  eroded  area,  Ae  =  eroded  cross-sectional  area,  yw  is  the  water  specific 
weight,  Nz  =  t/Tm  =  number  of  waves  at  mean  period  during  event  of  duration  t,  Tm  =  mean  wave 
period,  sm  =  HJLm  =  wave  steepness,  fls  significant  wave  height,  Lm  wavelength  based  on  mean  wave 
period,  and  9  is  the  seaward  structure  slope  from  horizontal. 

These  stability  equations  are  similar  to  those  proposed  by  Hudson  and  Van  der  Meer.  The  greatest 
difference  is  the  inclusion  of  a  maximum  momentum-flux-based  wave  force  derived  from  first 
principles.  The  maximum  wave  force  based  on  wave  momentum  flux  should  be  reasonably  well 
predicted,  even  in  shallow  water.  Previous  stability  relations  had  no  clear  dependence  on  water 
depth.  The  new  stability  relations  illustrate  the  well  known  fact  that  annor-  stone  stability  decreases 
for  increasing  wave  steepness  in  shallow  water  and  is  at  a  minimum  for  severe  plunging  breakers. 
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For  plunging  breakers,  stability  is  only  mildly  a  function  of  permeability  and  is  not  dependent  on 
wave  steepness.  For  surging  waves,  stability  is  more  strongly  related  to  permeability  and  wave 
steepness. 

These  equations  are  intended  for  preliminary  design.  It  is  recommended  to  utilize  physical  models, 
if  at  all  possible,  to  verify  stability  for  final  design.  Local  experience  is  valuable,  but  details  of  wave 
focusing  and  instability  at  structure  transitions  will  only  be  revealed  in  a  physical  model  study. 

EXAMPLE:  COMPUTE  STABLE  ARMOR  WEIGHT 

Find:  The  maximum  wave  momentum  flux  parameter  and  corresponding  stable  armor  weight  for  a 
conventional  nonovertopped  rubble-mound  breakwater  with  the  following  characteristics. 

Given: 

yw  =  64  pcf  -  water  specific  weight 

Sr  =2.6  -  stone  specific  gravity 

g  =  32.2  ft/s  -  acceleration  of  gravity 

Hs  =  7  ft  -  significant  wave  height  at  structure  toe 

Tm  =  12  sec  -  mean  wave  period 

h  =  14  ft  -  local  water  depth  at  structure  toe 

S  =  2  -  start  of  damage 

Nz  =  3,000  -  stonn  duration  =  6  hr  at  the  peak  of  the  storm 

P  =  0.4  -  structure  permeability  for  traditional  multilayer  cross  section 

cot  9  =  2  -  seaward  structure  slope  of  1V:2H 

Compute  the  maximum  wave  momentum  flux  at  the  structure  toe  using  Equation  4: 

/  \  2.026 

A  =0.639  —  =0.1569 

1 14  J 

V— 0.391 

4=0.180  —  =0.2360 

14  J 


4 


ERDC/CHL  CHETN-1 1 1-7 1 
December  2005 
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/max 


0.6172 


0 M? ) 


64(14)2(0. 1569) 


14 


-0.236 


32.2(12) 


7742.4 


Compute  local  steepness  of  mean  wavelength:  sm  =  HJLm  where  Lm  =  250  ft  is  the  local  mean 
wavelength  found  using  the  linear  wave  dispersion  relation  with  Tm  =  12  sec  and  h  =  14  ft. 
Therefore,  sm  =  7/250  =  0.0280. 

Compute  critical  wave  steepness  using  Equation  8:  smc  =  -0.0035(2)  +  0.028  =  0.021. 

In  this  case,  the  waves  are  plunging  on  the  structure  because  sm  >  smc.  Equations  6  and  9,  therefore, 
should  be  used  to  size  the  stone.  Equation  6  for  plunging  waves  on  the  structure,  yields: 

Nm  =  5 .0(2  /  3000° 5  )0-2  o.4018  yfl  =  3.093 


Then  the  stable  annor-stone  size  is  given  by  Equation  9  as 


D.  = 


1(0.6172) 


(2.6-1) 


1/2 


14 


3.093 


=  2.81  ft 


The  corresponding  armor  weight  is  W=  D„  Sr  yw  =  2.8 13(2. 6)64  =  3,698  lb. 

If  all  conditions  are  the  same  but  the  water  depth  is  28  ft,  then  the  calculations  would  proceed  as 
follows: 
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0.639 


2.026 


0.0385 
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0.180 


0.3095 


(■ ) 
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28 
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Compute  local  steepness  of  mean  wavelength:  sm  =  HJLm  where  Lm  =  346  ft  is  the  local  mean 
wavelength  found  using  the  linear  wave  dispersion  relation  with  Tm  =12  sec  and  h  =  28  ft. 
Therefore,  sm  =  7/346  =  0.0202. 

Compute  critical  wave  steepness  using  Equation  8:  smc  =  -0.0035(2)  +  0.028  =  0.021. 

In  this  case,  the  waves  are  surging  on  the  structure  because  sm  <  smc.  Equations  7  and  9,  therefore, 
should  be  used  to  size  the  stone.  Equation  7  for  surging  waves  on  the  structure,  yields: 

Nm  =  5.0(2/3000°'5)°'20.40'18  V2(.0202)"4/3  =5.202 


Then  the  stable  stone  size  is  given  by  Equation  9  as 


D 


1(0. 1873)  l‘  2  28 
(2.6  —  1)  5.202 


1.84  ft 


*3  T 

The  corresponding  armor  weight  is  given  by  W=  Dn  Sryw=  1 .84  (2.6)64  =  1,039  lb.  Comparingthe 
stone  weights,  the  shallow  water  depth  resulted  in  an  increase  in  the  stable  stone  weight  by  a  factor 
of  more  than  3.5.  The  majority  of  the  impact  of  shallower  depth  came  from  the  momentum  flux 
parameter  in  the  stability  number  calculation. 


POINTS  OF  CONTACT:  This  CHETN  was  developed  within  the  Risk  Analysis  of  Coastal 
Structures  work  unit  in  the  Navigation  Systems  R&D  Program.  The  program  is  administered  by 
Mr.  James  Clausner  of  the  Coastal  and  Hydraulics  Laboratory,  U.S.  Army  Engineer  Research  and 
Development  Center.  Questions  about  this  Technical  Note  can  be  addressed  to  Dr.  Jeffrey  A. 
Melby  (Jeflfev.A.Melbv@erdc.usace.armv.mil)  or  Dr.  Steven  Hughes  (Steven. A.Hughes@erdc 
.uscae.army.mil).  Questions  about  the  R&D  program  should  be  addressed  to  Mr.  Clausner 
(James.E.Clausner@erdc.usace.army.mil). 
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